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The Weber-Davis model of the solar wind is generalized to include the effects of latitude. ABSTRACT 

The principal assumptions of high electrical conductivity, rotational symmetry, the 

polytropic relation between pressure and density, and a flow-alined field in a system 

rotating with the sun, are retained. An approximate solution to the resulting equations 

for spherical boundary conditions at the base of the corona indicates a small component 

of latitudinal flow toward the solar poles at large distances from the sun as a result of 

latitudinal magnetic forces. 


The purpose of this paper is to demonstrate that even 
with spherically symmetric boundary conditions at the 
base of the corona, the magnetic field introduces a small 
latitude effect in the solar wind that should be 
measurably significant at large distances. As motivation 
for a more general development it is useful to consider 
the magnetic forces predicted by the spiral field 
approximation of Parker [1958]. The magnetic force 
per unit mass for this model is 

f V X BJX B/7 47rp y = -V [(B 0 2 /A'np a ) sin 2 0 ] (1) 

where B 0 and p Q are reference values for the magnetic 
field and density. The total force per unit mass is then 
given by the negative gradient of the sum of the thermal 
energy per unit mass, the gravitational energy per unit 
mass, and the magnetic energy per unit mass (eq. (1)). 
The thermal term is assumed to be falling adiabatically 
toward zero at large distances, and the gravitational term 
is decreasing toward zero; however, for spherical 
boundary conditions on B 0 and p Q , the magnetic term 
depends upon sin 2 9 and does not decrease explicitly 
with r. This suggests that changes in momentum at large 
distances may be dominated by the magnetic field. 
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The magnetic torque per unit mass for this model is 
then 

r= rX(dV/dt) = ~(B () 2 /4iTp 0 ) sin 6 cos 6 $ (2) 

which acts to drive the plasma toward the poles, 
producing a nonzero component of the 6 component of 
velocity Vq. For infinite electrical conductivity, the 
presence of a finite component of Vq implies a finite 
component of Bq through the process of convection. 
The assumed field is therefore inconsistent with respect 
to latitude and requires a more general treatment. It 
should be noted that the magnetic stress given by this 
model results because of the <p component of the 
magnetic field. Since this winding of the field is due to 
solar rotation, the effect demonstrated here may be 
expected to persist in a more general treatment. 

The approach to be described here is identical with 
that of Weber and Davis [1967] except that the 
equations are generalized to include the effects of 
latitude through the development of the normal 
component of momentum. We begin with the 
steady-state flow equations subject to the following 
assumptions: (1) neglect viscosity, (2) isotropic thermal 
conductivity, (3) infinite electrical conductivity, 

(4) equal electron and proton temperatures, and 

(5) energy equation may be replaced by the polytropic 
law. We further assume rotational symmetry and that 



the magnetic field and flow velocity are parallel in a 
reference system rotating with the sun. A flow-alined 
field may be generally represented by 

B =KpU (3) 

where U = V -coyr sin 9 $ is the flow velocity in the 
rotating system, V is the flow velocity in a fixed system 
of references is the mass density, and coy is the solar 
rotation rate. The condition V- B = 0, combined with 
the equation of continuity, requires k to be a streamline 
constant. The physical meaning of k is obtained by 
noting that the quantity k 2 pf 4rr must be unity at the 
Alfven point, or k 2 =4 where p^ is the density at 
the Alfven point. 

With rotational symmetry the </> component of the 
momentum equation may be directly integrated to 
obtain V ^ as a function of r, 9, and p, which when 
combined with the magnetic field as a function of 
velocity and density (eq. (3)), reduces the problem to 
one of two-dimensional flow in the rO plane. Given 
boundary conditions, the unknowns at a general point r, 
6 are p, V r , and Vg, which are determined by the 
equation of continuity, the parallel component of 
momentum (Bernoulli’s equation), and the normal 
component of momentum. The problem is complicated 
by the presence of three critical surfaces corresponding 
to the fast, slow, and Alfven modes of propagation in 
the medium. 

A particular streamline will be designated by 
specifying its coordinate 0 o at the base of the corona, 
which is constant along the streamline. There are six 
required boundary values that must be specified or 
determined as functions of 6 0 at the coronal base: 

T 0 (9 0 >• Pj e o>- v n/°o)’ v do(°o>- K ( d o>< and the total ^ 
component of angular momentum along the streamline 
L(9 0 ). The requirement that the velocity components be 
finite at the three critical surfaces supplies three of the 
required boundary conditions for Vg 0 , «, and L. 
Further, if BJ9 0 j is specified, then V ro may be 
determined through equation (3). Hence, we chose to 
specify T 0 (9 0 ), pJ9 0 ), and B o (0 o ), the magnitude of 
the total magnetic field at the coronal base. 

We- find that the best way to cope with the critical 
surfaces is to integrate the equations along a streamline 
until the critical point is encountered and then to vary a 
boundary condition until the streamline that passes 
through is determined. Because of the differential form 
of the equations, integration along a streamline is 
possible only if the previous neighboring streamline is 
known. Hence, to start the integration, the polar limit of 
the equations is required. However, the equations at 


the pole are not closed and an assumption is necessary to 
obtain the polar limit. A complete solution to the 
problem is possible through a process of iteration over the 
entire set of streamlines; however, techniques we have 
been able to devise thus far are too costly. 

As a first approximation to the solution, we chose to 
close the equations with the assumption 

v d/ v r = ~{ r W(r)/dr] !f(r)\ sin 9 cos 9 (4) 

which integrates to 

tan 9 0 =f(r) tan 9 (5) 

Equation (4) is the lowest order 9 dependence for which 
Vg is zero at both the equator and the pole. The radial 
function f(r) specifies the shape of the streamlines. 
Radial flow in the r9 plane is given by f(r) = 1 and 
df(r)/dr = 0. The function f(r) is determined by integra- 
tion of the limit of the normal component of momentum 
along the polar stream line. In the polar limit we have 

9 VqI'69 = -V r [r df(r)jdr] / f(r ) (9 = 0) (6) 

d9 0 /d9 =f(r) (7) 

The results of equations (4) and (5) and the polar 
integration of f(r) for two sets of constant temperature, 
density, and magnetic field at the coronal base are 
presented in tables 1 through 4 and figure 1. Figure 1 
shows the radial dependence of Vg along a streamline 
for the two sets of boundary conditions; tables 1 and 2 
tabulate the resulting values at 1 AU, and tables 3 and 4 
tabulate the corresponding critical surfaces. 

The function f(r) obtained at the pole is not overly 
sensitive to reasonable changes in the assumed model to 
close the equations for constant boundary conditions. 
The general shape of the curve results from the physics 
of the problem and is maintained. Since Vg for 0 o = 10° 
may be obtained directly from the polar solution 
through equation (6), this curve should be in close 
agreement with an exact solution. The curves for 
9 0 > 10°, however, are contingent upon how well the 
assumption contained in equation (4) approximates the 
normal component of momentum at lower latitudes. 

The results presented here indicate a modest 
latitudinal flow toward the poles due to magnetic 
stresses, which persist out to very large distances from 
the sun for spherical boundary conditions at the corona. 
Associated with this flow is a 9 component of the 
magnetic field. The radial flow parameters are not 
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Figure 1. Latitudinal flow velocity for two sets of constant boundary conditions at the 
coronal base, r CJ = I.OXiO 1 1 cm. Curves for 6 0 = 50° , 60°, 70° and 80° are nearly 
identical with the displayed curves for 40°, 30°, 20° and 10°, respectively, and are not 
shown. 
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Table 1 . Predicted values at 1 AU 

[Constant coronal boundary conditions: T Q = 2X10 6 °K ,N g = 3X10 6 particles/cm 3 , B 0 = 1 gauss] 


#0 > 
deg 

9, 

deg 

Velocity, km/sec 

Magnetic field, 10 5 

gauss 

Density , 
protons/cm 3 

Temperature, 

°K 

Vr 

Ve 

V 4> 

B r 

B e/ B r 

B (f>l B r 

10 

9.88 

403 

-0.28 

0.19 

3.57 

-7.06X1 O' 4 

-0.19 

13.6 

3.16X10 5 

20 

19.78 

403 

-0.54 

0.37 

3.57 

-1.33X10-* 

-0.38 

13.6 

3.16X10 5 

30 

29.70 

403 

-0.72 

0.53 

3.56 

-1.80X10 -3 

-0.55 

13.6 

3.16X10 5 

40 

39.65 

403 

-0.83 

0.68 

3.55 

-2.05X10 -3 

-0.71 

13.5 

3.16X10 5 

50 

49.65 

403 

-0.83 

0.81 

3.54 

-2 .06X10“ 3 

-0.85 

13.5 

3.16X10 5 

60 

59.69 

403 

-0.73 

0.91 

3.53 

-1.82X10 -3 

-0.96 

13.5 

3.15X10 5 

70 

69.77 

404 

-0.55 

0.98 

3.52 

-1.35X10-3 

-1.04 

13.5 

3.15X10 5 

80 

79.88 

404 

-0.29 

1.03 

3.51 

-7.22X1 0 -4 

-1.09 

13.5 

3.15X10 5 


Table 2. Predicted values at 1 AU 

[Constant coronal boundary conditions: T a = 2X10 6 °K, N 0 = 1X10 6 particles/ cm 3 , B 0 = 1.5 gauss) 


9 o' 

deg 

9, 

deg 

Velocity, km/sec 

Magnetic field, 10 

gauss 

Density, 

protons/cm 3 

Temperature, 

°K 

V r 

V B 


B r 

B e/ B r 

B <j>/ B r 

10 

9.48 

402 

-2.37 

0.60 

6.75 

-5.89X1 0~ 3 

-0.18 

5.85 

3.28X10 5 

20 

19.01 

402 

-4.50 

1.16 

6.70 

-1.12X icr 2 

-0.36 

5.83 

3.28X10 5 

30 

28.66 

403 

-6.15 

1.64 

6.64 

-1.53X10 -2 

-0.53 

5.80 

3.28X10 5 

40 

38.46 

403 

-7.13 

2.05 

6.56 

-I.77X10' 2 

-0.69 

5.75 

3.27X1 0 5 

50 

48.44 

403 

-7.27 

2.34 

6.47 

-1.80X10 -2 

-0.83 

5.71 

3.27X10 5 

60 

58.62 

404 

-6.51 

2.56 

6.39 

-1.61X10 -2 

-0.94 

5.66 

3.27X1 0 s 

70 

68.98 

404 

-4.81 

2.69 

6.32 

-1.19X 10~ 2 

-1.03 

5.62 

3.26X1 0 5 

80 

79.45 

404 

-2.64 

2.78 

6.27 

-6.54X1CT 3 

-1.09 

5.59 

3.26X1 0 s 





Table 3. Critical surfaces 

[Constant coronal boundary conditions: T 0 - 2X10 6 °K, 
N 0 = 3X10 6 particles/ cm 3 , B a = 1 gauss] 



Slow mode 

Alfven mode 

Fast mode 

0 o > 

AS, 

VS’ 

G4> 


r fi 

e f> 

deg 

solar 

radii 

deg 

solar 

radii 

deg 

solar 

radii 

deg 

10 

4.012 

9.997 

15.233 

9.979 

15.239 

9.979 

20 

4.011 

19.994 

15.226 

19.961 

15.236 

19.961 

30 

4.010 

29.992 

15.211 

29.947 

15.233 

29.947 

40 

4.008 

39.991 

15.194 

39.940 

15.230 

39.940 

50 

4.006 

49.991 

15.174 

49.940 

15.226 

49.940 

60 

4.004 

59.992 

15.159 

59.947 

15.224 

59.947 

70 

4.003 

69.994 

15.144 

69.961 

15.221 

69.961 

80 

4.002 

79.997 

15.136 

79.979 

15.220 

79.979 


Table 4. Critical surfaces 

[Constant coronal boundary conditions: T a = 2X10 6 °K, 
N 0 = IX 10 6 particles/cm 3 , B a = 1.5 gauss] 



Slow mode 

Alfven mode 

Fast mode 

®o> 


9 S ’ 

r A> 


r f’ 


deg 

solar 

radii 

deg 

solar 

radii 

deg 

solar 

radii 

deg 

10 

4.015 

9.997 

36.043 

9.895 

36.067 

9.895 

20 

4.013 

19.994 

35.931 

19.803 

36.027 

19.803 

30 

4.011 

29.992 

35.773 

29.736 

35.977 

29.734 

40 

4.008 

39.991 

35.578 

39.701 

35.914 

39.698 

50 

4.003 

49.991 

35.354 

49.703 

35.827 

49.698 

60 

4.001 

59.992 

35.178 

59.740 

35.781 

59.734 

70 

3.999 

69.994 

35.020 

69.805 

35.731 

69.803 

80 

3.997 

79.997 

34.931 

79.896 

35.703 



79.895 


significantly different from the Parker or Weber-Davis 
models for latitudinal flow velocities up to about 
10 km/sec. 
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